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INTRODUCTION 


The  rate  of  decomposition  of  an  ablative  material  when  heated  is  modeled  by  the  kinetic 
rate  equation.  If  it  is  assumed  that  the  material  dimensions  are  constant,  the  rate  equation  de¬ 
termines  the  density  of  the  remaining  char.  Both  the  rate  of  decomposition  and  char  density 
strongly  affect  the  thermal  performance  of  the  material.  In  order  to  predict  the  thermal  response, 
accurate  values  of  the  kinetic  parameters  over  the  entire  range  of  decomposition  are  required.  In 
the  case  of  an  ablative  material  exposed  to  a  solid  rocket  motor  exhaust,  the  heat  flux  may  vary 
widely  depending  upon  the  geometry  and/or  type  of  motor.  Therefore,  the  effect  of  the  heating 
rate  on  the  kinetic  parameters  must  be  known  as  well. 

The  purpose  of  this  study  was  to  determine  an  appropriate  model  for  the  rate  of  decom¬ 
position  of  ablative  materials.  Friedman’s  method1  using  multiple  heating  rates  was  chosen,  since 
ablative  materials  are  subjected  to  widely  varying  heating  rates.  Application  of  this  method  re¬ 
quired  calculating  an  average  activation  energy  for  the  entire  thermal  decomposition.  The 
decomposition  reaction  required  two  models,  one  for  the  initial  decomposition  and  another  for 
the  remainder.  For  these  two  regions,  separate  pre-exponential  factors  and  apparent  orders  of 
reaction  were  calculated  by  the  technique  developed  by  Friedman.1  For  comparison,  the  average 
activation  energy  for  each  material  determined  by  Friedman’s  method  was  compared  to  the  value 
obtained  by  the  method  of  Flynn  and  Wall.2 


BACKGROUND 

The  decomposition  kinetics  of  solid  materials  have  been  studied  by  many  investigators. 
As  a  result,  numerous  techniques  have  been  developed  to  extract  the  kinetic  parameters  from 
experimental  data. 

Freeman  and  Carroll3  developed  the  well-known  difference  method  and  applied  the 
technique  to  determine  the  kinetic  parameters  for  calcium  oxalate  monohydrate.  The  method 
was  later  revised  by  Anderson  and  Freeman4  and  applied  to  the  study  of  polystyrene  and 
polyethylene.  Mickelson  and  Einhom5  developed  the  ratio  method  to  analyze  thermogravimetric 
data  obtained  for  a  urethane  polymer.  Baer,  Hedges,  Seader,  Jayakar,  and  Wojcik6  heated  samples 
of  reinforced  polymers  at  heating  rates  up  to  4200°C/min.  The  data  were  correlated  by  a 
numerical  technique  developed  by  Bumingham  and  Seader.7  Friedman1  studied  the 
decomposition  of  a  fiberglass-phenolic  based  on  a  technique  developed  for  multiple  heating  rates. 
Similarly,  Flynn  and  Wall2  developed  a  method  for  determining  the  activation  energy  based  on 
data  taken  at  several  different  heating  rates. 

Baer,  et  al.6  discussed  the  fact  that  kinetic  parameters  obtained  by  methods  using  a  single 
thermogram  at  low  heating  rates  do  not  accurately  predict  kinetic  behavior  when  applied  to  the 
higher  heating  rates.  For  this  reason,  the  methods  of  Friedman1  and  Flynn  and  Wall2  were  con¬ 
sidered  in  this  work. 
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THEORY 


FRIEDMAN’S  METHOD 

Friedman’s  method  is  attractive  for  this  application  because  the  kinetic  properties  may  be 
calculated  based  on  data  taken  over  a  wide  range  of  heating  rates.  Further,  the  Arrhenius  equation 
is  combined  with  an  arbitrary  function  of  weight.  This  allows  more  flexibility,  since  no  prior 
knowledge  of  the  function  is  required.  This  method  does,  however,  require  measurement  of  the 
weight  loss  and  rate  of  weight  loss  as  a  function  of  temperature  at  several  different  heating  rates. 

The  general  form  of  the  rate  equation  proposed  by  Friedman  is 

- 1  /Wo  x  dW/dt  =  Af(W/Wo  )  exp  (-  E/RT)  ( 1 ) 


where 


W0 

dW/dt  = 

A 

E 

R 

T 

f(W/WQ)  = 


original  weight  of  material  (mg) 
rate  of  weight  loss  (mg/min) 
pre-exponential  factor  (min“ 1 ) 
activation  energy  (cal/gm-mole) 
gas  constant  (1.987  cal/gm-mole-°K) 
temperature  (°K) 
undefined  function  of  weight 


Taking  the  natural  logarithm  of  both  sides  of  Equation  (1 )  results  in 


ln[-l/Wo  x  dW/dt]  =  ln[Af(W/Wo)]  -  E/RT  (2) 

A  linear  equation  may  be  fit  to  ln(-  1/Wq  x  dW/dt]  as  a  function  of  1/T  at  constant  parametric 
values  of  W/WQ.  These  equations  will  have  slopes  of  -E/R.  Each  intercept  is  the  value  of 
ln[Af(W/Wo)]  at  the  parametric  value  of  W/WQ.  Then,  by  defining 

f(W/Wo)  =  [(W-  Wf)/WJn  (3) 


where 


n  =  order  of  reaction 
W  =  instantaneous  weight  of  material  (mg) 

Wf  =  final  weight  of  charred  material  (mg) 

and  multiplying  Equation  (3)  by  A  and  taking  the  natural  logarithm  results  in 

ln[Af(W/Wo)]  =  InA  +  nln((W- Wf)/Wo]  (4) 
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The  final  ratio,  Wf/WQ,  is  taken  from  the  original  thermograms.  Since  ln[Af(W/Wo)J  is  known 
for  various  W/WQ  ratios,  Equation  (4)  can  be  used  to  obtain  values  of  A  and  n. 


DISCUSSION  OF  FRIEDMAN  TECHNIQUE 

Friedman  used  this  technique  to  calculate  the  kinetic  properties  for  CTL91-LD  fiberglass- 
phenolic.  One  activation  energy  was  calculated  for  each  of  the  1 2  values  of  weight  loss  ranging 
from  0.675  to  0.95  (on  a  glass-free  basis).  The  average  activation  energy  was  calculated  from  these 
data.  By  eliminating  the  early  weight  loss  (a»4%)  and  dropping  the  data  points  above  W/Wo  = 
0.875  and  using  Wf/WQ  =  0.61 ,  a  linear  curve  was  fit  to  the  data.  Thus,  the  effective  range  covered 
by  the  curve  fit  was  approximately  0.65  <  W/Wo  <  0.85,  which  accounted  for  about  50  percent 
of  the  total  weight  loss.  This  resulted  in  a  rather  poor  fit  of  the  data  at  both  ends  of  the  weight 
loss  curve.  In  contrast,  the  requirement  for  the  present  application  was  to  obtain  a  kinetic 
expression  applicable  over  the  entire  range  of  weight  loss. 

Three  points  regarding  this  method  should  be  clarified.  First,  the  equation  of  f(W/Wo) 
may  take  a  variety  of  forms.  For  example,  Goldfarb,  McGuchan  and  Meeks8  selected 
f(W/Wo)  =  [(W  -  W{)/Wq  -  Wf)]n:  This  will  result  only  in  a  change  in  the  intercept  InA;  i.e.,  a 
change  in  the  apparent  pre-exponential  factor. 

Secondly,  the  kinetic  parameters  may  be  calculated  by  considering  either  the  total  weight 
or  only  the  resin  weight  of  the  sample.  Again  the  activation  energy  and  order  of  reaction  remain 
unchanged.  Only  the  intercept  InA  is  affected. 

Using  a  pre-exponential  factor  based  only  on  the  resin  weight  will  result  in  an  error,  if  used 
in  calculations  where  the  total  weight  is  being  considered.  Assuming  the  unknown  function  is 
[(W-  Wf)/WQ  ]",  the  two  pre-exponential  factors  are  related  by: 

A'  =  A[(Wq  -  Wg)/Wo)n-1  (5) 


where 


A'  =  pre-exponential  factor  (resin  weight  only)(min" 1 ) 

Wg  -  weight  of  inert  material  (mg) 

Finally,  changes  in  the  activation  energy  at  different  degrees  of  conversion  may  be  a 
result  of  real  changes  due  to  a  change  in  mechanism,  to  a  change  in  structure  of  the  resin,  or  to  a 
result  of  experimental  error.  If  these  changes  are  not  a  result  of  experimental  error,  then  using 
E  =  E(W/Wo)  would  be  more  realistic.  In  many  cases,  however,  separation  of  the  experimental 
error  from  real  changes  in  E  is  difficult. 


FLYNN  AND  WALL  METHOD 

Flynn  and  Wall2  developed  a  convenient  method  to  determine  the  activation  energy  from 
weight  loss  curves  measured  at  several  heating  rates.  The  following  relationship  is  used  to  calculate 
the  activation  energy. 
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E  w-  (R/C)dlog0/d(l/T) 


(6) 


where 


P  =  heating  rate  (°C/min) 

C  =  C(E/RT) 

Plotting  1/T  versus  log  p  at  several  weight  loss  ratios  results  in  a  series  of  straight  lines  with 
slope  Alog  PI A(l/T).  Using  the  slope  and  the  appropriate  value  of  C,  the  activation  energy  can  be 
calculated  by  Equation  (6).  Since  C  is  a  function  of  E/RT,  the  calculation  of  E  from  Equation  (6) 
is  an  iterative  process.  Elynn  and  Wall  constructed  a  table  of  values  for  C  over  the  range  from 
7  <  E/RT  <  60.  The  variation  of  C  over  this  range  is  approximately  ±3  percent.  This  method 
is  extremely  attractive,  since  it  involves  only  reading  the  temperature  at  a  constant  weight  loss 
from  a  series  of  thermograms  at  different  heating  rates. 


EXPERIMENTAL 


MATERIALS 

The  two  ablative  materials  studied  were  supplied  by  Haveg  Industries.  As  shown  in  Table  1 , 
these  materials  consisted  of  a  phenol-formaldehyde  resin  with  specified  amounts  of  glass,  asbestos, 
and/or  magnesium  silicate  added  as  filler. 

The  materials  were  converted  to  powder  form  by  machining  and  were  then  filtered  through 
a  No.  20  sieve.  They  were  stored  overnight  in  a  vacuum  dessicator  maintained  at  35°C  to  remove 
traces  of  water. 


Table  1.  Composition  of  Material  Tested 


Meforial  P nmnncifmn 


Contents 


Asbestos 
Glass  (Si02 )  and 
Magnesium  Silicate 
Total  Filler  Content 
Phenol-Formaldehyde 
Resin  (H4 IP) 

Total  Nonvolatiles 


H41NE 

H41D 

(%) 

(%) 

- 

52.0 

60.5 

_ 

60.5 

52.0 

39.5 

48.0 

60.5 

52.0 

4 


APPARATUS  AND  PROCEDURE 


A  Perkin-Elmer  TGS-2  Thermogravimetric  System  was  used,  with  temperature  control 
provided  by  a  Perkin-Elmer  System  4  Microprocessor  Controller.  The  sample  temperature  was 
measured  with  a  chromel-alumel  thermocouple,  which  was  calibrated  with  a  set  of  five  Curie 
standards  in  the  temperature  range  of  interest  at  each  heating  rate  used.9 

In  order  to  reduce  temperature  gradients  in  the  material  and  to  ensure  uniform  heating, 
small  weights  of  a  powdered  form  of  the  materials  were  used.  Samples  weighing  7.5±0.5  mg  were 
heated  from  40  to  950°C  using  heating  rates  of  10,  20,  40,  80,  100,  and  160°C/min.  Both  the 
percentage  of  initial  weight  and  the  rate  of  weight  loss  were  plotted  directly  as  a  function  of 
temperature.  The  samples  were  maintained  in  a  nitrogen  atmosphere  throughout  the  experiment. 
When  the  programmed  temperature  scan  reached  950°C,  the  purge  gas  was  automatically  switched 
to  oxygen  to  thermo-oxidatively  degrade  the  remaining  resin.  To  verify  the  initial  weight  fraction 
of  filler,  the  temperature  was  held  at  950°C  until  the  resin  had  completely  degraded. 


RESULTS 


The  original  thermograms  contained  the  temperature,  derivative  of  weight  loss,  and  the 
fraction  of  weight  remaining.  These  data  were  digitized  at  0.01  intervals  of  the  fraction  of  weight 
remaining.  The  experimental  temperatures  were  corrected  using  the  Curie  standard  temperature 
calibration  for  each  heating  rate.  The  thermograms  were  then  reproduced  from  this  data.  Figures  1 
through  6  and  7  through  12  show  the  fraction  of  weight  remaining  and  the  rate  of  weight  loss 
as  a  function  of  temperature  and  time  for  each  heating  rate  for  H41NE  and  H41D,  respectively. 
Comparisons  of  the  fraction  of  weight  remaining  as  a  function  of  temperature  at  all  six  heating 
rates  for  H41NE  and  H41D  are  shown  in  Figures  13  and  14.  Figures  15  and  16  depict  the  relative 
magnitudes  of  the  derivative  of  weight  loss  as  a  function  of  temperature  at  all  six  heating  rates 
for  H41NE  and  H41D,  respectively.  The  digitized  data  for  both  materials  are  listed  in  the 
appendix. 
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Figure  7.  160°C/min  Thermogram  for  H41D 
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Figure  8.  100°C/min  Thermogram  for  H41D 
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Figure  9.  80°C/min  Thermogram  for  H41D 
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Figure  10.  40°  C/min  Thermogram  for  H41D 
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Figure  11.  20°C/min  Thermogram  for  H41D 
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Figure  12.  10° C/min  Thermogram  for  H41D 


-1/Wq  x  DW/DT  -  1/MIN  x  103  -1/Wq  X  DW/DT  -  1/MIN  x  10* 


1.00 


HAVEG  41  NE 
WEIGHT  LOSS  FOR 
6  HEATING  RATES 


S  0.95 
O 

z 

z 

|  0.90 

ui 

K 


2  0.86 

ui 

i 

z 

o 

O  0.80  I- 

< 

K 


0.76 


J _ l _ i _ i _ i _ i _ i _ i - 1 

2  34  66789  10 

TEMPERATURE  -  °C  X10'1 


Figure  13.  Fraction  Weight  Remaining  for  All  Six  Heating  Rates  for  H41NE 


Figure  14.  Fraction  Weight  Remaining  for  All  Six  Heating  Rates  for  H41D 


A  plot  of  ln[-  1/WQ  x  dW/dt}  versus  I/T  for  H41NE  is  shown  in  Figure  17.  Figure  18  is  a 
a  single  plot  for  H41D.  The  slope  of  each  line  was  determined  from  a  least-squares  fit  of  the  data. 
Figure  19  shows  the  corresponding  activation  energy  and  intercept  ln[Af(W/W0)]  at  each  value  of 
weight  loss  from  0.84  <  W/WQ  <  0.98  for  H41NE.  Figure  20  depicts  the  same  information  for 
H41D  over  the  range  from  0.79  <  W/WQ  <  0.97.  The  range  of  each  data  point  is  the  range  of 
error  based  on  the  least-squares  fit  of  the  data.  Values  of  ln[Af(W/Wo)]  versus  ln[(W  -  Wf)/Wo  ] 
for  H41NE  are  shown  in  Figure  21  and  for  H41D  in  Figure  22.  These  figures  depict  the  separation 
of  the  reaction  into  the  two  regions  and  the  corresponding  least-squares  fit  over  each  region.  A 
pre-exponential  factor  and  order  of  reaction  were  determined  for  each  of  these  two  regions.  The 
average  activation  energies  determined  from  Figures  1 7  and  1 8  were  used  for  both  regions.  Using 
Flynn  and  Wall’s  method,  the  average  activation  energies  for  both  materials  were  calculated  based 
on  plots  of  log  0  versus  1/T.  The  results  are  shown  in  Figures  23  and  24  for  H41NE  and  H41D, 
respectively.  A  summary  of  the  results  of  the  calculations  for  both  materials  is  listed  in  Table  2. 

The  kinetic  parameters  calculated  by  the  modified  version  of  Friedman’s  method  were 
used  in  Equation  (1)  to  calculate  the  fraction  of  weight  remaining  versus  temperature.  Each  set 
of  parameters  was  applied  to  that  portion  of  the  weight  loss  curve  from  which  it  was  determined. 
A  comparison  of  the  results  of  these  calculations  and  the  experimental  data  for  10°C/min  and 
160°C/min  heating  rates  for  H41NE  and  H4ID  are  shown  in  Figures  25  and  26,  respectively. 
The  average  error,  standard  deviation  of  errors  and  the  95-percent  confidence  interval  were 
calculated  for  90  and  114  experimental  versus  calculated  points  for  H41NE  and  H41D, 
respectively.  These  results  are  presented  in  Table  3. 
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ure  17.  Plot  of  Slopes  Used  to  Determine  the  Activation  Energy  for  H41NE 


Figure  18.  Plot  of  Slopes  Used  to  Determine  the  Activation  Energy  for  H41D 
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Figure  19.  Activation  Energy  and  Intercept  as  a  Function 
of  Degree  of  Conversion  for  H41NE 
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Figure  20.  Activation  Energy  and  Intercept  as  a  Function 
of  Degree  of  Conversion  for  H41D 
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Figure  21.  Plot  to  Determine  the  Pre-exponential  Factor  and  Order  of 
Reaction  for  Two  Regions  of  Weight  Loss  for  H41NE 


Figure  22.  Plot  to  Determine  the  Pre-exponential  Factor  and  Order  of 
Reaction  for  Two  Regions  of  Weight  Loss  for  H41D 
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Figure  23.  Plot  of  Slopes  Used  to  Determine  the  Activation  Energy 
for  H41NE  by  the  Method  of  Flynn  and  Wall 
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Figure  24.  Plot  of  Slopes  Used  to  Determine  the  Activation  Energy 
for  H41D  by  the  Method  of  Flynn  and  Wall 


Table  2.  Summary  of  Calculations 

Range 


of 

Eg  v  K  (kcal/gm-mole) 

Material 

wf/w0 

W/W0 

Friedman 

Flynn  &  Wall 

A(min~ 1 ) 

n 

W/W0 

H41NE 

0.795 

0.98  -  0.84 

62.13 

62.15 

1.19  x  1031 

4.90  x  1020 

17.33 

6.30 

>  0.91 

<  0.91 

H41D 

0.760 

0.97  -  0.79 

74.97 

71.52 

2.71  x  103s 

3.87  x  1023 

19.46 

6.12 

>  0.89 

<  0.89 
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Figure  25.  Comparison  of  Calculated  Versus  Experimental  Weight  Loss  Errors 
for  H41NE  at  160°C/min  and  10°C/min 


Figure  26.  Comparison  of  Calculated  Versus  Experimental  Weight  Loss  Errors 
for  H41D  at  160°C/min  and  10°C/min  Heating  Rates 
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Table  3.  Statistical  Analysis  of  Errors  in  Computed 
Versus  Experimental  W/WQ 


Average 

Standard 

Error 

Deviation 

Confidence  Interval 

Material 

(%) 

(%) 

(95%) 

H41NE 

0.33 

0.58 

0.22  to  0.44 

H41D 

0.28 

0.84 

0.14  to  0.42 

DISCUSSION 


The  average  activation  energies  calculated  by  the  methods  of  Flynn  and  Wall  and  Friedman 
agree  within  0.03  and  4.6  percent  for  H41NE  and  H41D,  respectively.  There  was  less  scatter  of 
the  data  using  the  Flynn  and  Wall  method.  This  was  thought  to  be  due  primarily  to  the  errors  in 
measuring  the  derivatives  of  the  weight  loss  used  in  Friedman’s  method. 

The  kinetic  parameters  were  calculated  based  on  data  that  represented  68  percent  of  the 
total  weight  loss  for  H41D  and  75  percent  for  H41NE.  This  resulted  in  large  values  of  the  order 
of  reaction  and  pre-exponential  factor  for  the  first  region.  However,  Friedman’s  results  would  have 
been  similar  to  these  had  he  considered  the  same  range  of  decomposition. 

In  order  to  evaluate  the  effect  of  separating  the  reaction  into  two  parts,  the  thermograms 
were  calculated  for  H41NE  using  only  the  kinetic  parameters  for  W/Wq  <  0.91.  The  same 
calculations  were  made  for  H41D  for  W/WQ  <  0.89.  This  corresponded  approximately  to  the 
region  of  weight  loss  considered  by  Friedman.  As  shown  by  the  broken  lines  in  Figures  25  and  26, 
the  calculated  versus  the  experimental  thermograms  are  in  poor  agreement. 

By  separating  the  reaction  in  this  manner,  the  reaction  order  and  pre-exponential  factor 
become  empirical  parameters  that  provide  a  “best  fit”  of  the  data.  However,  this  method  yields 
an  extremely  accurate  reproduction  of  the  thermograms  over  a  wide  range  of  heating  rates.  This  is 
the  desired  result  for  kinetic  parameters  used  in  thermal  models. 
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Experimental  Data 


Experimental  Data  (Continued) 
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